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Abstract 

Observations are reported on isotactic polypropylene (i) in a series of tensile 
tests with a constant strain rate on specimens annealed for 24 h at various tem- 
peratures in the range from 110 to 150 °C and (ii) in two series of creep tests in 
the sub-yield region of deformation on samples not subjected to thermal treatment 
and on specimens annealed at 140 °C. A model is developed for the elastoplas- 
tic and nonlinear viscoelastic responses of semicrystalline polymers. A polymer 
is treated an equivalent transient network of macromolecules bridged by junctions 
(physical cross-links, entanglements and lamellar blocks). The network is assumed 
to be highly heterogeneous, and it is thought of as an ensemble of meso-regions 
with different activation energies for separation of strands from temporary nodes. 
The elastoplastic behavior is modelled as sliding of meso-domains with respect to 
each other driven by mechanical factors. The viscoelastic response is attributed to 
detachment of active strands from temporary junctions and attachment of dangling 
chains to the network. Constitutive equations for isothermal uniaxial deformation 
are derived by using the laws of thermodynamics. Adjustable parameters in the 
stress-strain relations are found by fitting the experimental data. 

1 Introduction 

This paper is concerned with the effect of annealing on the elastoplastic and viscoelas- 
tic responses of injection- molded isotactic polypropylene (iPP). This semi-crystalline 
polypropylene is chosen for the analysis because of its numerous applications in industry 
(oriented films for packaging, reinforcing fibres, nonwoven fabrics, blends with thermo- 
plastic elastomers, etc.). 

The effect of annealing at elevated temperatures on the morphology of polypropylene 
has been a focus of attention in past five years, see, e.g., |], §, ||, |], [5], §, [7J , to mention a few. 
The previous works concentrated on calorimetric studies of morphological transformations 
driven by primary and secondary crystallization. The influence of changes in the micro- 
structure of spherulites on the mechanical behavior of iPP was not studied in detail. 
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The nonlinear viscoelastic response of polypropylene was analyzed by Ward and Wolfe 
and Smart and Williams || several decades ago and, in the past decade, by Wortmann 
IJ, Ariyama [fT2]j , Ariyama et al. Dutta and Edward []3~4] , Read and 
|T6[ 1 and Sweeney et al. 



and Schulz 
Tomlins UTol 
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Tomlins and Read 

Yielding and viscoplasticity of iPP have been recently investigated by Aboulfaraj et 
TEfl , Kalay and Bevis [JJ, Coulon et al. [[B|, Seguela et al. PH , Staniek et al. pi| , 



5J, Seguela et al. PH , Staniek et al. 
Nitta and Takayanagi p2| , |23| and Labour et al. f24l , to mention a few. 

Isotactic polypropylene is a semi- crystalline polymer containing monoclinic a crystal- 
lites, hexagonal (3 structures, orthorhombic 7 polymorphs and "smectic" mesophase ||. 
At rapid cooling of the melt (at the stage of injection molding), a crystallites and smectic 
mesophase are mainly developed, whereas (3 and 7 polymorphs are observed as minority 
components [[l], [|. 

A unique feature of a spherulites in iPP is the lamellar cross-hatching: development 
of transverse lamellae oriented in the direction perpendicular to the direction of radial 
lamellae |5|, [|. The characteristic size of a spherulites in injection-molded specimens is 



19]. These spherulites consist of crystalline lamellae with 



estimated as 100 to 200 /im 
thickness of 10 to 20 nm |^ 

The amorphous phase is located (i) between spherulites, (ii) inside spherulites, in 
"liquid pockets" between lamellar stacks [^5| , and (iii) between lamellae in lamellar stacks. 
It consists of (i) relatively mobile chains between spherulites, in liquid pockets and between 
radial lamellae inside lamellar stacks, and (ii) severely restricted chains in the regions 



bounded by radial and tangential lamellae (rigid amorphous fraction [25]). 

Stretching of iPP specimens results in inter-lamellar separation, rotation and twist of 



lamellae, fine and coarse slip of lamellar blocks and their fragmentation [IS], 20|, chain slip 
through the crystals, sliding, pull-out and breakage of tie chains |22|, ^3j, and activation 



of rigid amorphous fraction. At large strains, these morphological transformations lead 



to cavitation, formation of fibrills and stress- induced crystallization of iPP |26 



To develop stress-strain relations, we apply a method of "homogenization of micro- 
structure," according to which a sophisticated morphology of isotactic polypropylene is 
modelled by an equivalent phase whose deformation captures essential features of the 
mechanical response. An amorphous phase is chosen as the equivalent phase for the 
following reasons: 

1. The viscoelastic response of isotactic polypropylene is conventionally associated with 
rearrangement of chains in amorphous regions [TP] . 



Sliding of tie chains along and their detachment from lamellae play the key role in 
the time-dependent response of iPP [E], |2l 



3. The viscoplastic flow in semi-crystalline polymers is assumed to be "initiated in the 
amorphous phase before transitioning into the crystalline phase" P7[ . 



4. 



The time-dependent behavior of polypropylene is conventionally modelled within 

u 



the concept of a network of macromolecules 

Dynamic mechanical analysis reveals that the loss tangent of iPP demonstrates two 
pronounced maxima being plotted versus temperature [^, p9| . The first maximum 
transition in the interval between T = —20 and T = 10 °C) is associated with the glass 
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transition in a mobile part of the amorphous phase, whereas the other maximum (a- 
transition in the interval between T = 70 and T = 110 °C) is attributed to the glass 
transition in the remaining part of the amorphous phase, the rigid amorphous fraction 
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At room temperature (i.e., above the glass transition temperature for the mobile 
amorphous phase) iPP is treated a transient network of macromolecules [3C] bridged 



by junctions (physical cross-links, entanglements and lamellar blocks). The network is 
assumed to be highly heterogeneous (the inhomogeneity is attributed to interactions be- 
tween spherulites and amorphous regions, as well as to local density fluctuations in the 
amorphous phase), and it is thought of as an ensemble of meso-regions (MRs) with dif- 
ferent activation energies for separation of strands from temporary nodes. 

Two types of MRs are distinguished: (i) active domains, where strands separate from 
junctions as they are thermally agitated (mobile part of the amorphous phase), and (ii) 
passive domains, where detachment of chains from junctions is prevented by surrounding 
lamellae. 

Stretching of a specimen induces 

1. slippage of meso-domains with respect to each other (which models the elastoplastic 
behavior of iPP), 

2. separation of active strands from temporary junctions and attachment of dangling 
chains to the network in active meso-regions (which reflects the viscoelastic re- 
sponse), 

3. an increase in the content of active MRs driven by release of the rigid amorphous 
fraction due to lamellar fragmentation (which is associated with the nonlinearity of 
the viscoelastic behavior in the sub- yield region of deformation). 

The objective of this study is three-fold: 

1. To report experimental data in tensile tests with a constant strain rate on specimens 
annealed at various temperatures in the interval between 110 and 150 °C. 

2. To present obsevations in creep tests on non-annealed samples and on specimens 
annealed at 140 °C in the sub-yield region. 

3. To derive constitutive equations for the time-dependent behavior of a semicrystalline 
polymer and to find adjustable parameters in the stress-strain relations by fitting 
observations. 

The goal of this paper is to shed some light on correlations between morphological trans- 
formations in iPP at annealing and changes in stress-strain diagrams and creep curves 
measured in the sub-yield and post-yield regions of deformation. 

The exposition is organized as follows. The specimens and the experimental procedure 
are described in Section 2. Section 3 deals with kinetic equations for rearrangement of 
strands in active meso-domains. Sliding of MRs with respect to each other is described 
in Section 4. Stress-strain relations are derived in Section 5 by using the laws of ther- 
modynamics. Adjustable parameters in the constitutive equations are found in Section 6 
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by fitting observations. A brief discussion of our findings is presented in Section 7. Some 
concluding remarks are formulated in Section 8. 

2 Experimental procedure 

Isotactic polypropylene (Novolen 1100L) was supplied by BASF (Targor). ASTM dumb- 
bell specimens were injection molded with length 148 mm, width 10 mm and height 3.8 
mm. Mechanical tests were performed on non-annealed specimens, as well as on sam- 
ples annealed in an oven for 24 h at the temperatures 110, 120, 130, 140 and 150 °C. 
After annealing, the specimens were slowly cooled by air. To minimize the effect of phys- 
ical aging on the time-dependent response, tests were carried out a week after thermal 
pre-treatment. 

Uniaxial tensile tests were performed at room temperature on a testing machine 
Instron-5568 equipped with electro-mechanical sensors for the control of longitudinal 
strains in the active zone of samples (with the distance 50 mm between clips). The tensile 
force was measured by a standard load cell. The longitudinal stress, a, was determined 
as the ratio of the axial force to the cross-sectional area of stress-free specimens. 

In the first series of tests, non-annealed specimens and specimens annealed at various 
temperatures, T, were loaded with the cross-head speed 5 mm/min (that corresponded 
to the Hencky strain rate e# = 1.1 ■ 10 -3 s _1 ) up to the maximal strain e max = 0.06. 
According to the chosen strain rate ensured nearly isothermal test conditions. The 
engineering stress, a, was plotted versus the longitudinal strain, e, in Figure 1. 

A series of 6 creep tests was performed on non-annealed specimens at the longitudinal 
stresses a? = 10.0 MPa, a° = 15.0 MPa, a° = 17.5 MPa, a° 4 = 20.0 MPa, a° = 22.5 
MPa and erg = 25.0 MPa. Each creep test was carried out on a new sample. In the mth 
test (m = 1, . . . , 6), a specimen was loaded with the cross-head speed 5 mm/min up to 
the engineering stress cr^ that was preserved constant during the creep test, t c = 20 min. 
The longitudinal strain, e, was plotted versus the logarithm (log = log 10 ) of time t (the 
instant t = corresponds to the beginning of a creep test) in Figure 2. 

To evaluate the effect of annealing on the nonlinear viscoelastic behavior of iPP, a 
series of 5 creep tests was performed on specimens annealed at 140 °C at the longitudinal 
stresses a? = 10.0 MPa, <j% = 15.0 MPa, cx° = 20.0 MPa, al = 25.0 MPa, and trg = 30.0 
MPa. The strain, e, was plotted versus the logarithm of time, t, in Figure 3. 

Figures 2 and 3 demonstrate that the rate of increase in strain, e, with time, t, is 
relatively low at small stresses, and it noticeably grows with the longitudinal stress, a. Our 
aim is to derive constitutive equations for the elastoplastic and viscoelastic responses of a 
semi-crystalline polymer and to find adjustable parameters in the stress-strain relations 
by fitting the experimental data plotted in Figures 1 to 3. 

3 Rearrangement of strands in active meso-regions 

A semicrystalline polymer is treated as a transient network of macromolecules bridged by 
temporary nodes. The network is modelled as an ensemble of meso-regions with various 
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potential energies for detachment of strands from their junctions. Two types of meso- 
domains are distinguished: 

1. passive, where all nodes are thought of as permanent, 

2. active, where active strands (whose ends are connected to contiguous junctions) 
separate from the nodes at random times when these strands are thermally agitated. 

An active chain whose end detaches from a junction is transformed into a dangling chain. 
A dangling chain returns into the active state when its free end captures a nearby junction 
at a random instant. 

Denote by X a the number (per unit mass) of active strands in active MRs, and by X p 
the number (per unit mass) of strands connected to the network in passive MRs. Under a 
time-dependent loading program, some lamellae (restricting mobility of chains in passive 
MRs) break, which results in an increase in the number of strands to be rearranged. The 
quantities X a and X p are treated as functions of time, t, that obey the conservation law 

X Bt (t)+X v (t)=X, (1) 

where X is the average number of active strands per unit mass of a polymer (which is 
assumed to be time-independent). 

Rearrangement of strands in active MRs is thought of as a thermally activated process. 
The rate of detachment of active strands from their junctions in a MR with potential 



energy uj is given by the Eyring equation [32 



r = raeXP (~^ 

where &b is Boltzmann's constant, T is the absolute temperature, and the pre-factor 
T a is independent of energy uj and temperature T. Confining ourselves to isothermal 
deformations at a reference temperature To and introducing the dimensionless activation 
energy uj = uj/ (k^T ), we arrive at the formula 

r = r a ex P (-^). (2) 

An ensemble of active MRs with various potential energies is described by the distribution 
function p(t, uj) that equals the ratio of the number, N a (t, uj), of active meso-domains with 
energy uj at instant t to the total number of active MRs, 

N (t uj) r°° 

PM = -^jir, *.(*)=/ N*(t,u)du, (3) 
X a {t) Jo 

and by the concentration of active MRs 

«.(*> = M 

In what follows, constitutive equations will be derived for an arbitrary distribution func- 
tion p(t,uj). To fit experimental data, a random energy model is applied with 



p(t,uj) =p (t)exp 



2S 2 (t) 



(u)>0), p(t,u) = (u<0), (5) 
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where f2 is the average activation energy in an ensemble of active meso-domains, E is the 
standard deviation of potential energies for separation of strands, and po(t) is determined 
by the condition 

poo 

/ p(t,uj)duj = 1. (6) 

J 

An ensemble of active meso-domains is characterized by the function n a (t, r, uj) that equals 
the number (per unit mass) of active strands at time t belonging to active MRs with 
potential energy uj that have last been rearranged before instant r e [0, £]. In particular, 
n a (t,t,uj) is the number (per unit mass) of active strands in active MRs with potential 
energy uj at time t > 0, 

n & (t,t,u) = N a (t,uj). (7) 

The amount ip{r,uj)dr ) where 



On 



(8) 



t=T 



equals the number (per unit mass) of dangling strands in active MRs with potential energy 
uj that merge with the network within the interval [r, r + dr] , and the quantity 

— {t,r,uj)dr 

is the number of these strands that have not detached from temporary junctions during 
the interval [r,t\. The number (per unit mass) of strands in active MRs that separate 
(for the first time) from the network within the interval [t, t + dt] reads 

a -(t,0,u)dt. 



dt 

The number (per unit mass) of strands in active MRs that merged with the network 
during the interval [r, r + dr] and, afterwards, separated from the network within the 
interval [t,t + dt] is given by 

~dt^^ T '^ dtdT - 

The rate of detachment, T, equals the ratio of the number of active strands that separate 
from the network per unit time to the current number of active strands. Applying this 
definition to active strands that merged with the network during the interval [r, r + dr] 
and separated from temporary junctions within the interval [t, t + dt], we find that 

^(t,r,uj) = -T(uj)^(t,r,uj). (9) 

Changes in the function n a (t, 0, a;) are governed by two processes: 

1. detachment of active strands from temporary nodes, 

2. mechanically-induced activation of passive MRs. 
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The kinetic equation for this function reads 

dn 



dN a 



(t, 0, u>) = -T(u)n a (t, 0, u) + —^(t, u). 



The solution of Eq. (10) with initial condition (7), where we set t — 0, is given by 
n a (t, 0, cj) = iV a (0, uj) exp 



-r(w)*l + / 

It follows from Eqs. (8) and (9) that 

dn 



'* ON, 



dt 



(r, uj) exp 



-r(u)(t-r) 



dr. 



(t,T,u) = <p(T,u)exp[-r(u)(t-T) . 
Equation (7) implies that 

rt Qfi 

N a (t,u) = n a (t,0,u) +/ -—^(t,T,u)dr. 

Jo or 

Differentiating Eq. (13) with respect to time and using Eq. (8), we find that 



<p(t, u) + —(t, 0, u) + ^ — (t, r, u;)rfr = w ). 

This equality together with Eqs. (9), (10) and (13) results in 

<p(t,w) = T(u)N a (t,u). 
This expression together with Eq. (12) yields 

dn. 



dr 



■(t,T,u) = T(uj)N Si (t,uj)exp -r(u)(t - t) 



(10) 



fin 



(12) 



(13) 



(14) 



Rearrangement of strands in active MRs is described by Eqs. (2), (11) and (14). Sep- 
aration of active strands from their junctions and detachment of dangling chains to the 
network reflect the viscoelastic response of a semi-crystalline polymer. 



4 Sliding of meso-regions 

It is assumed that meso-domains are not rigidly connected, but can slide with respect 
to each other under straining. Sliding of meso-domains is treated as a rate-independent 
process and is associated with the elastoplastic behavior of a semi-crystalline polymer. 
We suppose that an increase in strain, e, by an increment, de, causes growth of the 
elastoplastic strain, e p , by an increment, de p , that is proportional to de, 

de p = ^/de. 

The coefficient of proportionality, depends, in general, on the macro-strain, e, the 
macro-stress, a, and the elastoplastic strain, e p . We presume that \1> is a function of the 
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elastic strain, e e , which is defined as the difference between the macro-strain, e, and the 
elastoplastic strain, e p , 

e e (t) = e(t)-e p (t). (15) 
This results in the kinetic equation 

^(t) = *(e(t)-e p (t))|(t), e p (0)=0. (16) 

It is natural to suppose that the function ^(e e ) vanishes at e e = (no elastoplastic strains 
are observed at small macro-strains), monotonically increases with the elastic strain, and 
reaches some limiting value b e (0, 1) at large values of e e (which corresponds to a steady 
regime of plastic flow). To minimize the number of adjustable parameters in constitutive 
equations, an exponential dependence is adopted, 



exp(— ae e 



(17) 



where the positive coefficients a and b are found by matching observations. It is worth 
noting that Eqs. (16) and (17) differ from conventional flow rules in elastoplasticity, where 
the elastoplastic strain is assumed to be proportional to stress, a. These equations are, 
however, in good agreement with the conclusion by Men and Strobl |33| that "tensile 
deformations of semi- crystalline polymers . . . are strain-controlled." 



5 Constitutive equations 

An active strand is modelled as a linear elastic medium with the mechanical energy 

w(t) = i/ie 2 (t), (18) 

where [i is the average rigidity per strand and e is the strain from the stress-free state to 
the deformed state of the strand. 

For strands belonging to passive meso-domains, the strain e coincides with the elastic 
strain e e . Multiplying the strain energy per strand, Eq. (18), by the number of strands 
in passive MRs, we find the mechanical energy of meso-domains where rearrangement of 
chains is prevented by surrounding lamellae, 

W p (t) = ^X p (t)e 2 c (t). (19) 



With reference to the conventional theory of temporary networks p0| , stresses in dangling 
strands are assumed to totally relax before these strands merge with the network. This 
implies that the reference (stress-free) state of a strand that is attached to the network 
at time r coincides with the deformed state of the network at that instant. For active 
strands that have not been rearranged until time t, the strain e(t) coincides with e e (t), 
whereas for active strands that have last been merged with the network at time r G [0, t], 
the strain e(t, r) is given by 

e(t,r) = e c (t) - e c (r). 
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Summing the mechanical energies of active strands belonging to active MRs with various 
potential energies, u, that were rearranged at various instants, r e [0, t], we find the strain 
energy of active meso-domains, 



W a (t) = -n J Q \n a (t,0,u)e 2 c (t) + J o -^(t,T,u)[e c (t)-e c (r) 

The mechanical energy per unit mass of a polymer reads 

W(t) = W^(t)+W p (t). 
It follows from this equality and Eqs. (15), (19) and (20) that 

W(t) = ^{x p (t)( e (t)-e p (t)) 2 + ^ OO [n a (t,0,^)( e (t)-e p (t)) 



dr \duj. 



(20) 



+ 



J^(t,r,co)((e(t)-e p (t)) - (e(r)-e p (r))) 2 dr 



duj 



(21) 



Our aim now is to calculate the derivative of the function W(t) with respect to time t. 
Differentiation of Eq. (21) results in 



— {t)= A( t )[-( t )--£(t)}-Mt), 



(22) 



where 



A(t) = n{Xp(t) (e(t) - ep(f)) + J o °° [n a (f , 0, u) (e(t) - e p (t)) 



+ 



j^(t,r,u;)((e(t)-e p (t)) - (e(r) - e p (r)))dr 



duo 



(23) 



Mt) = -^{^ E W(^)-^W) 2 + / o 00 [^(^o,c)( e (t)- ep (t)) 



+ 



f*^(t,T,u,)((e(t)-e p (t)) - (e(r) -e p (r))) 2 dr]du]. (24) 



It follows from Eqs. (1), (3), (13) and (23) that 



A{t) = ix 



X(e(t)-e p 



duo 



o dr 



(t,T,u)(e(r) - e p (r))o?T 



(25) 



Substituting Eqs. (9) and (10) into Eq. (24) and using Eqs. (1) and (3), we obtain 



1 r°° 

A 1 (t) = -fi / T(uo)doo 

2 Jo 



n 



Xt,0,u)(e(t)-e p (t)) : 



+ J Q -^(t,r,cv)((e(t)-e p (t)) - (e(r)-e p (r))) 2 dr 
For isothermal uniaxial deformation, the Clausius-Duhem inequality reads 

Qit)- 
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(26) 



dW , . 1 . s de . s 

._ w + _ ffW _ w >o, 



where p is mass density, and Q is internal dissipation per unit mass. Substitution of Eqs. 
(16) and (22) into this equation results in 



a(t)-pA(t)(l-*(e(t)-e p (t))) 



j t (t)+pA 1 (t)>0. 



(27) 



It follows from Eq. (26) that the function A\(i) is non-negative for an arbitrary program 
of loading. This means that the dissipation inequality (27) is satisfied, provided that the 
expression in the square brackets vanishes. This condition together with Eq. (25) results 
in the stress-strain relation 



a(t) 



E 



1 - ¥(e(f) - ep(f)) 
* dn 



(e(f) - e P 

J du J Q ^r(*> r > w )( e ( r ) - e p( T )) dr 



where E = ppX is an analog of Young's modulus. Substitution of Eqs. (3), (4) and (14) 
into this equality implies that 



a(t) = E[l-V(e(t)-e p (t))]\(e(t)-e p (t)) - K a (f) f°° T{uj)dw 

J 



x 



J Q exp(-r(w)(t-r))(e(r) - e p (r))p(r, u)dr 



(28) 



Given functions p(t,u) and K a (t), Eqs. (2), (16), (17) and (28) determine the time- 
dependent behavior of a semicrystalline polymer at isothermal uniaxial deformation. 

To approximate the experimental data reported in Section 2, we concentrate on tensile 
tests with constant strain rates and on creep tests. 

Confining ourselves to "rapid" tensile tests, when the effect of the material's viscosity 
on the mechanical response may be disregarded, we neglect the integral term in Eq. (28) 
and arrive at the constitutive equation 



a(t) = E 



1 _ *( e (f) _ ep ( f ))] ( e (f) _ €py 

Equations (16), (17) and (29) are determined by 3 material parameters: 

1. the elastic modulus E, 

2. the constant a that characterizes the rate of elastoplastic strain, 

3. the constant b that describes a developed plastic flow. 

It is natural to assume that in a standard creep test with a longitudinal stress a , 



(29) 



a(t) 



0, t<0, 
a , t > 0, 



the quantities « a and p are functions of the stress intensity a . It follows from Eqs. (15) 
and (28) that the elastic strain, e e , reads 



e e (f) 



1 " *(e e (*)) 



+ 



kJct ) / Z(t,uj)p(a°,uj)duj, 
Jo 



(30) 
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Z(t,uj) =T(uj) / exp — T(uj)(t — r) e c (r)dr. 



where 



Equation (31) implies that the function Z(t,u>) satisfies the differential equation 
dZ 



dt 



(t,u) = T(u) e e (t) - Z(t,u) 



Z(0,u) = 0. 



(31) 



(32) 



After solving Eqs. (30) and (32), the longitudinal strain, e, is found from Eqs. (15) and 

(16) which can be presented in the form 

=[l- *(e c (t))p^(t), 6(0) = eV), te(O) = e° e (a°), (33) 

where the initial conditions, e°(cr°) and €®(<r°), are determined by integration of Eqs. (16), 

(17) and (29) from a = to a = a . 

Given a stress a , Eqs. (2), (5), (30), (32) and (33) are characterized by 4 adjustable 
parameters: 

1. the average potential energy for separation of active strands Q, 

2. the standard deviation of potential energies of active MRs E, 

3. the concentration of active meso-domains « a , 

4. the attempt rate for detachment of strands from temporary junctions T a . 

It follows from Eqs. (2), (5), (30) and (32) that the quantities Q and T a are mutually 
dependent: an increase in fl results in an increase in T a . To reduce the number of 
material constants to be found by matching observations in creep tests, we set T a = 1 
s. Our purpose now is to find adjustable parameters E, a, b, K a , Q and E by fitting 
experimental data depicted in Figures 1 to 3. 



6 Fitting of observations 

We begin with the approximation of the stress-strain curves presented in Figure 1. Under 
uniaxial tension with the cross-head speed 5 mm/min, the strain e max = 0.06 is reached 
within 69 s. According to Figures 2 and 3, changes in strain induced by rearrangement of 
active strands during this period are insignificant at stresses up to 20 MPa, whereas the 
duration of stretching at higher stresses does not exceeds 30 s, which causes rather small 
growth of strains. Based on these observations, we treat the deformation process as rapid 
and apply Eqs. (16), (17) and (29) to fit experimental data. 

For any temperature of annealing, T, the stress-strain curve is matched independently. 
To find the constants E, a and b, we fix the intervals [0,a max ] and [0,6 max ], where the 
"best-fit" parameters a and b are assumed to be located, and divide these intervals into J 
subintervals by the points = iAa and bj = jAb — 1, . . . , J) with Aa = a max / J and 
Ab = b max /J. For any pair, {a^, &.,•}, we integrate the governing equations numerically 
(with the step Ae = 5.0 ■ 10~ 5 ) by the Runge-Kutta method. The elastic modulus 
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E = E(i,j) is found by the least-squares algorithm from the condition of minimum of the 
function 

where the sum is calculated over all experimental points, e m , depicted in Figure 1, a cxp 
is the longitudinal stress measured in a tensile test, and <7 num is given by Eq. (29). The 
"best-fit" parameters a and b minimize K on the set {cii,bj = 1,..., J)}. After 

determining their values, and bj, this procedure is repeated twice for the new intervals 
[aj_i,aj+i] and [bj-i, bj+i] to ensure an acceptable accuracy of fitting. 

The "best-fit" parameters E, e* = a -1 and b are plotted in Figures 4 to 6 as functions 
of the annealing temperature T. The experimental data are approximated by the linear 
functions 

E = E + EiT, e* = e* + e*iT, b = b + b ± T, (34) 

where the coefficients E m , e* m and b m (m = 0, 1) are determined by the least-squares 
technique. 

We proceed with fitting the creep curves for non-annealed specimens depicted in Figure 
2. For any stress, a , the quantities e° and e° are found by integration of Eqs. (16), (17) 
and (29) with the material constants found in the approximation of the stress-strain 
curve plotted in Figure 1. To determine fi(cr°), E(cr ) and «(<7°), the following algorithm 
is employed. We fix the intervals [0, Q max ], [0, £ max ] and [0, « max ], where the "best-fit" 
parameters Q, E and k are assumed to be located, and divide these intervals into J 
subintervals by the points Qj = iAQ, £j = jAE and K k = kAn k = 1, . . . , J) with 
AQ = fimax/J, AS = S max /J and An = K m3X /J. For any pair, {Qj,Ej}, the constant 
p° = p°(i,j) is found from Eq. (6), where the integral is evaluated by Simpson's method 
with 200 points and the step Au; = 0.15. For any triple, Ej, n k } : Eqs. (2), (5), 
(30), (32) and (33) are integrated numerically (with the time step At = 0.1) by the 
Runge-Kutta method. The "best-fit" parameters Q, E and k minimize the function 

fc) ^ ' j^cxp(^m) ^num(^m) 5 
tm 

where the sum is calculated over all experimental points, t m , presented in Figure 2, e exp 
is the strain measured in the creep test, and e num is given by Eq. (33). After determining 
the "best-fit" values, Jlj, Ej and K/~, this procedure is repeated for the new intervals 
[Ej_!,Ej + i] and [/Cfc-i, /c^+i], to provide an acceptable accuracy of fitting. 
The adjustable parameters Q, E and k are plotted versus the engineering stress, a, in 
Figures 7 to 9. The experimental data are approximated by the linear functions 

Q = Q + Qxcr, E = E + ExO", k = k + K±a, (35) 

where the coefficients Q m , E m and K m (m = 0, 1) are found by the least-squares technique. 

Finally, we approximate the experimental data in creep tests on specimens annealed 
for 24 h at T = 140 °C presented in Figure 3. To determine adjustable parameters in 
the constitutive equations, the same procedure of fitting is used as for the observations 
depicted in Figure 2. The "best-fit" quantities f2, E and k are plotted versus the stress a 
in Figures 7 to 9 together with their approximations by Eqs. (35). 



12 



7 Discussion 



Figure 1 demonstrates fair agreement between the observations in tensile tests with a 
constant strain rate and the results of numerical simulation. 

Figure 4 shows that the elastic modulus E monotonically increases with annealing 
temperature T. In the region of temperatures from room temperature to 130 °C, this 
increase is rather small, but the rate of growth in E(T) substantially increases in the 
interval of temperatures between 130 and 150 °C. This conclusion is in agreement with 
observations in calorimetric tests by other researchers, which reveal that the melting 
temperature, T m , monotonically grows with crystallization temperature, T, and the slope 
of the curve T m (T) noticeably increases at the temperature T = 130 °C. Results of three 
DSC (differential scanning calorimetry) studies on iPP are presented in Figure 10. In this 
figure, the dependence of the melting peak, T m , on the crystallization temperature, T, is 
aproximated by the linear function 

T m = c + ciT. (36) 

The coefficients c m (m = 0, 1) in Eq. (36) are found by using the least-squares algorithm. 

According to the Gibbs-Thomson theory (see, e.g., || 0]), an increase in the equi- 
librium melting temperature, T m , is tantamount to an increase in the average lamellar 
thickness. Comparing Figures 4 and 10, we draw a conclusion that the growth of the elas- 
tic modulus of iPP with annealing temperature may be attributed to lamellar thickening. 

Figures 5 and 6 reveal that the plastic strain, e*, and the rate of developed plastic 
flow, b, decrease in the interval from 110 to 130 °C and increase at higher temperatures. 
Two possible explanations may be provided for the decrease in e*: 

1. According to Labour et al. ||24|| , annealing of isotactic polypropylene in the range 
of temperatures between 110 and 130 °C results in an increase in the concentration 
of ductile f3 spherulites that enhance plastic flow. This phenomenon is explained 
by the fact that "in the hexagonal phase, the chains are loosely packed, which 
suggests that chains are mobile" f34|. At annealing above 130 °C, development of 



(3 spherulites is thermo dynamically unfavorable (Al-Raheil et al. |2j reported that 
"no (3 spherulites were found above the crystallization temperature 132 °C"), which 
implies an increase in e* and a corresponding decrease in elastoplastic deformations. 



2. Men and Strobl |33| recently suggested that development of subsidiary lamellae 
at annealing is a two-step process consisting of formation of lamellar blocks at 
the initial stage and their aggregation into a "blocky substructure" at the final 
stage. According to this scenario, enchancement of plastic flow at low-temperature 
annealing is attributed to sliding of isolated lamellar blocks formed at secondary 
crystallization (the duration of thermal pre-treatment at temperatures in the range 
of 110 to 130 °C is assumed to be insufficient for the development of a blocky 
structure), whereas its decay at annealing above 130 °C is ascribed to aggregation 
of these blocks into stacks of relatively rigid lamellae. 

Both explanations are rather far from being exhausted. The first is based on the as- 
sumption that j3 spherulites grow at low-temperature annealing, which was questioned in 
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several works, see, e.g., []5| and the references therein. The other explanation presumes a 
new mechanism for development of subsidiary lamellae, which has not yet been confirmed 
experimentally. It is worth also noting that these models do not establish links between 
changes in the plastic strain, e*, and appropriate alternations in the rate of plastic flow, 
b. 

A decrease in the rate of steady plastic flow, 6, with temperature of annealing be- 
low 130 °C is attributed to lamellar thickening, which results in slowing down of the 
developed plastic flow (which is associated with lamellar fragmentation and motion of 
isolated lamellar blocks). An increase in b at higher temperatures is explained by the fact 
that at high-temperature annealing two processes occur simultaneously: (i) thickening of 
dominant lamellae and growth of subsidiary lamellae and (ii) annihilation of transverse 
lamellae. Gu et al. [f?j reported that annealing of iPP in the high-temperature region leads 
to "the complete absence of the cross-hatching." This implies a substantial reduction in 
the rigidity of spherulites, and, as a consequence, an enhancement of the developed plastic 
flow. 

To assess the level of elastoplastic strains at stretching of iPP specimens annealed 
at different temperatures, we integrate Eqs. (16) and (17) with adjustable parameters 
determined by matching experimental data in tensile tests. The results of numerical 
simulation are presented in Figure 11. This figure reveals that the elastoplastic strain, e p , 
is practically independent of annealing temperature (although the stress-strain diagrams 
depicted in Figure 1 are strongly affected by thermal pre-treatment). 

Figure 7 demonstrates that the average activation energy for separation of strands from 
temporary nodes, Q, decreases with stress, a, for non-annealed samples and increases 
for specimens annealed at 140 °C. This qualitative difference in the dependence fi(cr) 
for annealed and non-annealed iPP may be associated with transformation of smectic 
mesophase into lamellar blocks during thermal treatment. In a stress-free non-annealed 
specimen, the concentration of flocks of smectic mesophase ("arrays of chains with a 
better order in longitudinal that in transverse chain direction" is rather large. These 
ordered arrays of macromolecules slow down rearrangement of strands in active MRs, 
which implies that at small stresses, the average activation energy of samples not subjected 
to thermal treatment is rather large. According to the Men-Strobl model , annealing 
of specimens results in transformation of smectic mesophase into lamellar blocks and 
subsidiary lamellae. This transformation of "clusters of slow junctions" in amorphous 
regions into the crystalline phase enhances the rearrangement process, which implies that 
at small stresses, Q for non-annealed specimens substantially exceeds that for annealed 
samples. 

With an increase in stress, clusters of smectic mesophase in amorphous regions are 
disintegrated, and these "ordered" chains are transformed into ordinary ones. This pro- 
cess enchances detachment of strands from temporary junctions, which is reflected by a 
decrease in the average potential energy of active meso-domains. 

On the contrary, the growth of longitudinal stress results in fragmentation of "weak" 
lamellar blocks formed from smectic mesophase at annealing. Pieces of broken lamellar 
blocks distributed in amorphous regions serve as extra physical cross-links with high 
activation energy for separation of strands. As a result, stretching of annealed specimens 
increases the potential energy for detachment of chains and slows down the rearrangement 
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process. 

This scenario is confirmed by experimental data for the standard deviation of potential 
energies of active meso-domains depicted in Figure 8. This figure reveals that £ decreases 
with stress for non-annealed specimens (which is attributed to destruction of clusters of 
smectic mesophase in amorphous regions) and increases with stress for annealed samples 
(which is ascribed to the growth of concentration of extra physical cross-links induced by 
breakage of "weak" lamellar blocks). 

For the quasi-Gaussian distribution function (5), Q and S may be treated as an appar- 
ent average potential energy for detachment of active strainds, and an apparent standard 
deviation of potential energies for separation of strands from the network. The average 
activation energy, Q Q , and the standard deviation of activation energies, S , are given by 



Jo 

We determine the dimensionless quantities Q and S according to these formulas and 
calculate the ratio 

that characterizes the width of the quasi-Gaussian distribution. The parameter £ is plotted 
versus stress, a, in Figure 11. The experimental data are approximated by the linear 
function 

$ = £o + 6T, (37) 

where the coefficients £ m (m = 0, 1) are found by the least-squares technique. Figure 11 
demonstrates that the width of the distribution of active MRs, Eq. (5), is practically not 
affected by thermal treatment, and it weakly increases with longitudinal stress. The rate 
of increase in £ with stress appears to be independent of annealing temperature. 

Figure 9 reveals qualitatively different effects of the longitudinal stress on the concen- 
tration of active meso-regions, k, for non-annealed specimens and for specimens subjected 
to thermal pre-treatment. For non-annealed specimens, k remains constant up to er* 18 
MPa and linearly decreases with stress at higher longitudinal stresses. For specimens, 
annealed at 140 °C, k slightly increases with stress in the entire interval of deformations 
under consideration. 

At relatively small stresses, o < <j*, the fraction of active meso-domains in non- 
annealed specimens noticeably exceeds that for annealed samples. This conclusion is 
explained by the fact that annealing of iPP results in secondary crystallization of a part 
of the amorphous phase, which implies that the content of active MRs is reduced. An in- 
crease in k with stress for annealed specimens is also quite natural, because it is associated 
with mechanically-induced fragmentation of weak lamellae and release of the amorphous 
phase, whose rearrangement was prevented by surrounding lamellae in a stress-free spec- 
imen. 

To explain a pronounced decrease in k for non-annealed specimens at stresses exceeding 
a*, we should recall that the fitting procedure presumed the quantities Q, S and k in 
Eqs. (5) and (30) to be uniquelly determined by stress, a, and to be independent of 
elastoplastic strain, e p . To assess, whether this hypothesis is adequate, we integrate 
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numerically Eqs. (2), (5), (30), (32) and (33) with the adjustable parameters found by 
matching observations and calculate the elastoplastic strain, e p , as a function of time, t. 
The results of numerical simulation are depicted in Figures 13 and 14. 

Figure 14 shows that for all stresses (except for the highest stress, a = 30 MPa), the 
elastoplastic strain, e p , slightly increases with time in a fashion typical of primary creep. 
Curve 5 in this figure demonstrates a transition from the primary creep to the secondary 
creep (a linear increase in e p with time) at the final stage of the creep test. 

In contrast with results presented in Figure 14, Figure 13 reveals that e p increases 
with time following the primary-creep mode only at relatively small stresses (a < a*). At 
higher stresses, the primary creep is transformed into the secondary creep (curves 4 and 
5), and, finally, into the ternary creep (curve 6). 

Based on this observation, one can conclude that below the critical stress, a*, the 
concentration of active MRs, k, is independent of the elastoplastic strain, e p . This im- 
plies that in the interval of stresses [0,0"*], k remains constant (in agreement with the 
experimental data plotted in Figure 9). The independence of k of longitudinal stress, a, 
in this interval may be ascribed to the fact that at relatively small deformations smec- 
tic mesophase is mainly rearranged, while dominant lamellae (that restrict mobility of 
amorphous meso-domains) remain undisturbed. 

Above the critical stress, a*, i.e., in the region of secondary and ternary creep flows, 
the parameter k becomes a function of two arguments: the longitudinal stress, a, and the 
elastoplastic strain, e p . This implies that the data for n depicted in Figure 9 at a > a* 
should be treated as "average" (over the creep curves) quantities. This means that the 
noticeable decrease in k (curve lb in Figure 9) should be attributed to the influence of 
secondary and ternary creep flows on rearrangement of strands in active MRs. 

The following scenario may be proposed for a decrease in k observed in creep tests 
at relatively large stresses. A pronounced growth of elastoplastic strain e p (curves 4 to 6 
in Figure 13) is associated with fragmentation of dominant lamellae and plastic flow of 
lamellar blocks. In contrast with annealed specimens, where primary lamellae are broken 
into small pieces that serve as extra physical cross-links in active MRs (curves 2 in Fig- 
ures 7 and 8), dominant lamellae in specimens not subjected to thermal treatment are 
assumed to be disintegrated into relatively large blocks whose average length is compa- 
rable with the characteristic size of active MRs. Plastic flow of these blocks results in 
confining of some amorphous regions, where rearrangement of strands becomes prevented 
by surrounding immobile lamellae and moving lamellar blocks. This confinement of active 
MRs is temporary: when a lamellar block moves away from a "trapped" meso-region, it is 
released, and the rearrangement process proceeds. However, the higher the plastic strain, 
e p , is, the large is the number of "moving lamellar blocks," and, as a consequence, the 
smaller is the concentration of active MRs. 

According to this picture, an increase in the content of active meso-domains with stress 
for annealed samples is attributed to breakage of lamellae into small blocks, whereas a de- 
crease in k with stress for non-annealed specimens is ascribed to fragmentation of lamellae 
into large blocks, whose size is comparable with that of active MRs. This difference in the 
processes of lamellar disintegration for annealed and non-annealed samples is in agree- 
ment with the micro-mechanisms of lamellar growth at secondary crystallization recently 



proposed by Hikosaka et al. [34]. According to it, lamellar thickening is accompanying by 
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transformation of folded chain crystals in the lamellar cores into extended chain crystals. 
Packing of extended chains in the central part of a growing lamella becomes substantially 
looser compared to the initial one, which implies that dominant lamellae after thickening 
at secondary crystallization are fragmented into noticeably smaller pieces than the same 
lamellae before annealing. 

8 Concluding remarks 

A model has been developed for the elastoplastic and viscoelastic responses of semicrys- 
talline polymers at isothermal loading. A complicated micro-structure of a semi-crystalline 
polymer is replaced by an equivalent transient network of macromolecules bridged by 
junctions (physical cross-links, entanglements and crystalline lamellae). The network is 
thought of as an ensemble of meso-regions with various potential energies for separation 
of strands from temporary nodes. 

The viscoelastic response of a semicrystalline polymer is attributed to (i) detachment of 
active strands from temporary nodes in active meso-domains and (ii) merging of dangling 
strands with the network. Rearrangement of strands is treated as a thermally activated 
process, whose rate is determined by the Eyring equation (2). 

The elastoplastic behavior is ascribed to slippage of junctions with respect to their 
positions in the bulk material. Kinetic equation (16) is proposed for the rate of sliding, 
where the plastic flow is governed by the elastic strain (not stress, as in conventional 
theories of plasticity). 

Stress-strain relation (28) has been developed for isothermal uniaxial deformation by 
using the laws of thermodynamics. Adjustable parameters in the constitutive equations 
are found by fitting experimental data in tensile tests with a constant strain rate and in 
creep tests. 

A series of tensile tests has been performed on isotactic polypropylene at room tem- 
perature. The mechanical experiments are carried out on injection-molded specimens not 
subjected to thermal treatment and on samples annealed for 24 h at the temperatures 
110, 120, 130, 140 and 150 °C. 

Two series of creep tests have been performed on non-annealed specimens and on 
specimens annealed at 140 °C in the interval of stresses from 10 to 30 MPa. Fair agreement 
is demonstrated between the experimental data and the results of numerical simulation. 

The following conclusions are drawn: 

1. The flow rule (15), (16) and (29) with 3 adjustable parameters correctly describes 
observations in "rapid" tensile tests, when the effect of material viscosity is negligi- 
ble. The elastic modulus, E, increases with annealing time (which is attributed to 
lamellar thickening), whereas the plastic strain, e*, and the rate of developed plastic 
flow, b, demonstrate a more sophisticated behavior: they linearly decrease with an- 
nealing temperature, T, below T c = 130 °C and increase above T c . The temperature 
T c roughly coincides with the temperature, at which the rate of lamellar thickening 
noticeably increases. 
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2. The governing equations (2), (5), (30), (32) and (33) with 3 material constants 
adequately describe experimental data in creep tests under the assumption that 
the quantities Q, E and k are stress-dependent. The average activation energy for 
separation of chains from temporary nodes, Q, increases with stress for annealed 
specimens and decreases for samples not subjected to thermal treatment. This dif- 
ference is attributed to different mechanisms of lamellar fragmentation into blocks. 
The width of the quasi-Gaussian distribution (5) is weakly affected by longitudinal 
stress and annealing temperature. 

3. The concentration of active meso-domains, k, grows with stress for annealed spec- 
imens. This parameter is stress-independent for non-annealed samples below the 
threshold stress, a*, and decreases with stress at a > c». The latter observation 
is ascribed to transition from the primary creep of non-annealed specimens to the 
secondary and ternary creep flows, when rearrangement of strands in active MRs 
becomes strongly affected by elastoplastic strains. 
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List of figures 



Figure 1: The stress a MPa versus strain e in tensile tests. Circles: experimental data. 
Curve 1: a non-annealed specimen; curves 2 to 6: specimens annealed at a temperature 
T °C. Curve 2: T = 110; curve 3: T = 120; curve 4: T = 130; curve 5: T = 140; curve 6: 
T = 150. Solid lines: results of numerical simulation 

Figure 2: The strain e versus time t s in a tensile creep test with a stress o MPa. 
Circles: experimental data for non-annealed specimens. Solid lines: results of numerical 
simulation. Curve 1: a — 10.0; curve 2: o = 15.0; curve 3: o = 17.5; curve 4: a = 20.0; 
curve 5: a = 22.5; curve 6: a = 25.0 

Figure 3: The strain e versus time t s in a tensile creep test with a stress a MPa. Circles: 
experimental data for specimens annealed at T = 140 °C. Solid lines: results of numerical 
simulation. Curve 1: a = 10.0; curve 2: a = 15.0; curve 3: a = 20.0; curve 4: a = 25.0; 
curve 5: a = 30.0 

Figure 4: The elastic modulus E GPa versus annealing temperature T °C. Circles: 
treatment of observations. Solid lines: approximation of the experimental data by Eq. 
(34). Curve 1: E = 2.0242, E x = 0.0010; curve 2: E = 0.4946, Ex = 0.0129 

Figure 5: The plastic strain e* versus annealing temperature T °C. Circles: treatment 
of observations. Solid lines: approximation of the experimental data by Eq. (34). Curve 
1: e* = 3.53 • 10~ 2 , e* x = -1.08 • 10" 4 ; curve 2: e* = 8.68 • 10" 3 , e*i = 1.04 • 10" 4 

Figure 6: The rate of developed plastic flow b versus annealing temperature T °C. Circles: 
treatment of observations. Solid lines: approximation of the experimental data by Eq. 
(34). Curve 1: b = 0.9614, b x = -0.0018; curve 2: b = 0.4657, b x = 0.0021 

Figure 7: The average potential energy for detachment of strands Q versus stress a 
MPa. Symbols: treatment of observations. Unfilled circles: non-annealed specimens; filled 
circles: specimens annealed at T = 140 °C. Solid lines: approximation of the experimental 
data by Eq. (35). Curve 1: fi = 8.5079, fii = -0.1756; curve 2: fi = 4.3232, fi x = 
0.0950 

Figure 8: The standard deviation of potential energies for separation of strands £ ver- 
sus stress a MPa. Symbols: treatment of observations. Unfilled circles: non-annealed 
specimens; filled circles: specimens annealed at T = 140 °C. Solid lines: approximation 
of the experimental data by Eq. (35). Curve 1: S = 4.4272, S x = —0.0811; curve 2: 
S = 0.8672, Si = 0.1347 

Figure 9: The concentration of active MRs k versus stress a MPa. Symbols: treatment of 
observations. Unfilled circles: non-annealed specimens; filled circles: specimens annealed 
at T = 140 °C. Solid lines: approximation of the experimental data by Eq. (35). Curve 
la: k = 0.48, K X = 0.0; curve lb: k = 1.1052, k x = -0.0367; curve 2: k = 0.2371, 
Kl = 0.0057 

Figure 10: The melting peak T m °C versus crystallization temperature T °C. Symbols: 
treatment of observations for isotactic polypropylene. Unfilled circles: iPP crystallized for 
22 h at temperature T and cooled to room temperature 0; filled circles: iPP crystallized 
for an unspecified time at temperature T 0; asterisks: iPP crystallized for 100 min at 
temperature T 0. Solid lines: approximation of the experimental data by Eq. (36). 
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Vertical lines indicate transition temperatures 

Figure 11: The elastoplastic strain e p versus strain e in tensile tests on specimens an- 
nealed at various temperatures T. Solid lines: results of numerical simulation 

Figure 12: The ratio £ versus stress a MPa. Symbols: treatment of observations. Unfilled 
circles: non-annealed specimens; filled circles: specimens annealed at T = 140 °C. Solid 
lines: approximation of the experimental data by Eq. (37). Curve 1: £o = 0.4643, 
fi = 0.0014; curve 2: £ = 0.4369, ft = 0.0011 

Figure 13: The elastoplastic strain e p versus time t s in a tensile creep test with a stress 
a MPa. Solid lines: results of numerical simulation for non-annealed specimens. Curve 
1: a = 10.0; curve 2: a = 15.0; curve 3: a = 17.5; curve 4: a = 20.0; curve 5: a = 22.5; 
curve 6: a — 25.0 

Figure 14: The elastoplastic strain e p versus time t s in a tensile creep test with a stress 
a MPa. Solid lines: results of numerical simulation for specimens annealed at 140 °C. 
Curve 1: a — 10.0; curve 2: a = 15.0; curve 3: a = 20.0; curve 4: a = 25.0; curve 5: 
a = 30.0 
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